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ABSTRACT 

Structure formation provides a strong test of any cosmic acceleration model because a suc- 
cessful dark energy model must not inhibit or overpredict the development of observed large-scale 
structures. Traditional approaches to studies of structure formation in the presence of dark energy 
or a modified gravity implement a modified Press-Schechter formalism, which relates the linear 
overdensities to the abundance of dark matter haloes at the same time. We critically examine the 
universality of the Press-Schechter formalism for different cosmologies, and show that the halo 
abundance is best correlated with spherical linear overdensity at 94% of collapse (or observation) 
time. We then extend this argument to ellipsoidal collapse (which decreases the fractional time of 
best correlation for small haloes), and show that our results agree with deviations from modified 
Press-Schechter formalism seen in simulated mass functions. This provides a novel universal pre- 
scription to measure linear density evolution, based on current and future observations of cluster 
(or dark matter) halo mass function. In particular, even observations of cluster abundance in a 
single epoch will constrain the entire history of linear growth of cosmological of perturbations. 

Subject headings: cosmology: large-scale structure of universe 


1. Introduction 

An important part of the effort to explain cos- 
mic acceleration and the cosmological constant 
problem is testing proposed models in the context 
of what are, at this stage, better-established phys- 
ical pictures. Structure formation could prove to 
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be an incredibly useful phenomenological method 
for distinguishing models of cosmic acceleration. 

It is currently believed that large-scale struc- 
ture formation has its seeds in small quantum 


fluctuations in the early universe (e.g., Mukhanov 


et al. (1992)). The current model for structure 


formation is elegant in its fundamental simplicity. 
Random inhomogeneities, artifacts of cosmic in- 
flation, create a runaway effect in which overdense 
regions attract more matter, thus becoming more 
dense and leading to galaxies, stars, and planets. 
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Better understanding of this process is inde- 
pendently an intriguing enterprize in the field of 
cosmology. However, in this work we focus on 
the relationship between the cosmic acceleration 
and structure formation. More specifically, differ- 
ent cosmological pictures (cosmologies with differ- 
ing causes of acceleration, such as a cosmologi- 
cal constant, dark energy, or modifications of Ein- 
stein gravity) might have expansion histories that 
are similar to one another but leave different im- 
prints on large-scale structures, and in particular 
on galaxy clusters. Therefore, structure forma- 
tion provides a unique testing ground for models 


of cosmic acceleration (e.g., 

Islrak et al. 

(20061; 

Acquaviva et al. 

to 

O 

O 

OO 

Vikhlinin et al. ( 

2009 

))■ 


Here we critically examine some of the assump- 
tions in this program, and develop a framework to 
enhance the accuracy of this kind of work. 


The first step in this direction is to revisit how 


the Press-Schechter formalism (Press & Schechter 


1974) (PSF) is used to predict the cluster mass 


function. Press & Schechter (1974) have argued 


that the number density of dark matter haloes (or 
galaxy clusters) of mass M is given by: 


dn(M, z) p m {z)d lncr 1 (M,z) 

~^M~ = f[a(M,z)\— — , 

(1) 

where cr 2 (M, z) is the variance of linear overden- 
sity in spherical regions of mass M at redshift z, 
while p m {z) is the mean matter density of the uni- 
verse. For random gaussian initial conditions, f[cr\ 
is given by: 


fpsW } 



2 a 2 


(2) 


where 6 SC (~ 1.68 for most ACDM cosmologies) is 
the spherical collapse threshold for linear overden- 
sities (Gunn & Gott 1972). 

While the PSF successfully predicts the broad 
features of the simulated cluster mass functions, 
it proves too simplistic for detailed model com- 
parisons required for precision cosmology. Conse- 


quently, several authors including Sheth & Tormen 

(1999 

), Jenkins et al. (2001), 

Evrard et al. (2002), 

Warren et al. 

(2006), and Tinker et al. 

(2008 

have 


refined the function f(cr) to better match the mass 
functions in N-body simulations. We shall refer to 
these approaches as modified PSF. For example, 


Warren et al. (2006) and Tinker et al. (2008) pro- 


pose a fitting formula of the form: 


/(o') = A 


(!)' 
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( 3 ) 


where ( A,a,b,c ) = (0.186,1.47,2.57,1.19) give a 
good fit to simulated haloes of overdensity A = 
200, in a concordance ACDM cosmology at z = 0 
(Tinker et al. 2008) (see Fig. [l]). While most of 
this work is based on fitting formulae to simulated 
mass functions, Sheth et al. (2001) argue that an 
approximate implementation of ellipsoidal collapse 
can account for most of the deviations from the 
PSF. 



Fig. 1. - A comparison of Press-Schechter prediction for 
the function /(cr) (dotted line; Eq. [ 2 ]), with a parameterized 
fit to the numerical simulations (solid line; Eq. [3|. 


However, a more pressing question for cosmo- 
logical applications is whether the function /(cr) is 
universal, or rather can vary for different cosmolo- 
gies or cosmic acceleration models. In other words, 
could the same modified PSF be used to accurately 
predict halo abundance in cosmologies with differ- 
ent cosmological parameters? While earlier stud- 
ies failed to find such dependence, Tinker et al. 
(2008) first noticed a systematic evolution of /(cr) 
with redshift, implying a breakdown of universal- 


(2007); 


the 10 — 30% level (also see 

Reed et al. 

Crocce et al. 

2010); Bhattacharya et al. 

. Courtin et al. 

(2010) note that universal- 


ity is limited by the nonlinearity of structure for- 
mation, and the cluster mass function shows some 
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redshift dependence at higher redshifts that can be 
(partially) understood in the context of spherical 
collapse. However, the spherical collapse model 
falls short of explaining the 10 — 30% deviations 
from universality in all but the most massive clus- 
ters (see Fig. [7|. 

In this work, we contribute to the effort to bet- 
ter understand the role and limits of universal- 
ity in the cluster mass function by introducing a 
new parameter that appears to be universal across 
cosmological modelsj^] In particular, the modi- 
fied PSF relies on cr(M), the root-mean-square 
of linear density fluctuations at the time of ob- 
servation, when in reality, observed clusters are 
very non-linear objects with overdensities exceed- 
ing 200. We thus seek to find a universal time 
in the past when we could make a connection be- 
tween the nonlinear structures that we observe in 
the present and the linear structures that existed 
at that time, since all structures go through a lin- 
ear phase. Our basic strategy is to find the time in 
the past at which the linear density of the struc- 
tures that collapse today show minimum disper- 
sion as we vary cosmologies. This is illustrated in 
Fig. (2). 

In § [2j we present the complete nonlinear differ- 
ential equation that governs the growth of matter 
perturbations in spherical overdense regions in the 
presence of a cosmological constant. 

In § [3j we describe a numerical method which 
we developed in order to solve both the nonlin- 
ear and linear structure formation equations in 
the presence of a cosmological constant, or dark 
energy with a constant equation of state. 

In § |4j we discuss the implications of our nu- 
merical study, i.e. we find a universal fraction of 
collapse time, where linear density is the same in- 
dependent of cosmology. 

In § [5] we reexamine the idea of universality of 
mass functions in light of the results of the previ- 
ous section, including the effect of ellipsoidal col- 
lapse on the formalism, and propose a new mass 
function for general dark energy models. 


1 Here we define cosmological models to mean different val- 
ues of the cosmological constant density a t the time 
of observation, or different values of the dark energy equa- 
tion of state parameter, w. For now, we do not consider 
quintessence models where w is dynamical in nature, or 
possible modifications of Einstein gravity. Future work will 
extend our study to such possibilities. 


Finally, in §[bj we conclude with an overview of 
our results and a discussion of future prospects. 


2. A &: Non-linear Structure Formation 


In Appendix [Aj we derive the differential equa- 
tion that governs the growth of linear matter per- 
turbations. We used a familiar fluid dynamics pic- 
ture to do so. Here we derive the full non-linear 
equation for spherical overdensities using only cos- 
mological considerations. It should be noted that 
this particular form of the non-linear equation is 
only strictly valid for ACDM cosmologies, where 
the inside of a spherical top-hat overdensity can be 
considered as a separate universe. More complex 
models such as dark energy models with different 
values of w require additional considerations which 
will be discussed below. 

We consider a physical picture in which struc- 
ture formation arises due to a uniformly positive 
spherical perturbation away from an average mat- 
ter density, i.e. a top-hat matter overdensity. This 
scenario is similar to considering two cosmologies 
with two distinct scale factors: one for the outer 
universe and another for the overdense region. Of 
course, we are interested in a scenario where a 
dark energy component similar to a cosmological 
constant is at play, so we will assume the presence 
of one as part of our base model. 

For the external universe, we write the Fried- 
mann equation with zero curvature: 
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87 tGjv 
3 


(p + Pa) — H~. 


( 4 ) 


We can assume ~p oc a -3 for ordinary matter, while 
p\ = const, denotes the cosmological constant 
density. We also note that the value of the cosmo- 
logical constant will be the same inside the over- 
dense region and the background. 


In a general scenario with dynamical dark en- 
ergy models we cannot assume that curvature, of- 
ten denoted by k , will be a constant inside the 
overdense region due to the presence of pressure 
gradients. Therefore, Wang & Steinhardt (1998) 
point out that we are compelled to use the time- 
time component of Einstein’s equations, as these 
do not explicitly involve the curvature term. How- 
ever, in the presence of a cosmological constant, or 
w = — 1, we can ignore these considerations and 
begin with the first order Friedmann equation. We 
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Current Proposal 


Standard Universal Mass Functions 




Fig. 2. — A schematic picture of our proposal for a universal mass function: While standard universal mass functions relate 
the observed cluster abundance to linear growth at the same time, we propose that cluster abundance is best predicted by linear 
growth at an earlier time. The top charts show the abundance for clusters of different mass (nominally, masses of 10 13 , 10 14 
and 10 15 /i — 1 Mq, from top to bottom) as a function of cosmic time, while the bottom charts show the linear growth. As we 
argue in § [5] the abundance of less massive clusters depends on the linear growth at an earlier time, which explains why the 
arrows point to different stages in the history of linear growth. 


then calculate k, which can be seen as an integra- 
tion constant that arises in going from second to 
first order Friedmann formulations, using initial 
conditions. 

The scale factor in the overdense region is gov- 
erned by: 


8ttGn 


(, Pi + Pa)- 


( 5 ) 


Again, we define pi = p( 1 + 5). A little bit of al- 
gebra yields the following full differential equation 
for 6: 


[p(l + <5)] 3 


8t tG n 2 HS 

— V- o 


1 

+ tr- 


6 2 


3 1 + S 9(1 + S) 2 


+k = 0 


( 6 ) 

It is important to reiterate the importance of 
having access to both linear and nonlinear solu- 
tions. As noted by Pace et al. (2010) amongst 
others, although initially both the complete solu- 
tion and its linear approximation will track, even- 
tually the nonlinear solution will grow much faster 
relative to the scale factor. 


Following Liddle & Lyth (2000), we can show 
that knowing nonlinear theory is necessary. A crit- 
ical point in the evolution of a structure’s collapse 
is the turnaround event in which universal expan- 
sion’s dominance over the perturbation is eclipsed 
by gravitational collapse. In other words, at the 
turnaround point, a potential structure has de- 
tached from background expansion, but complete 
gravitationally-bound structure formation has not 
yet begun. This might be thought as the true birth 
of a structure within the void. Knowing the non- 
linear solution allows us to find out the value of the 
scale factor and the overdensity at this so-called 
“turn-around point.” 


Numerically, at the point of complete collapse 
the nonlinear solution will “blow up” and ap- 
proach infinity (see Fig. [3|. The linear density at 
this point, S sc is the quantity that enters the origi- 
nal Press-Schechter formalism (see §[!]) and is used 
as proxy between linear and non-linear structures. 
Physically, we do not expect an actual singularity. 
This “blow up” point is considered to be the be- 
ginning of virialization, a process whereby energy 
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in the bulk infall of matter is redistributed into 
random motion of dark matter particles, leading 
to a system in virial equilibrium, where kinetic en- 
ergy T and the potential energy U are related by 
the virial theorem: 


T, 


vir =2 & dU/dR) v 


( 7 ) 


(see Maor & Lahav ( 2005 ) ) . Another point of view 
on the spherical collapse picture, including a dis- 
cussion of virialization, can be found in |Somerville| 
& Primack (1999). However, for the purpose of 


this study, more details on this process are not 
necessary. 

We can generalize to cases where w ^ — 1, 
i.e. non-cosmological constant scenarios. |Abramo| 
et al. (2007) provides a comprehensive derivation 


of the full non-linear equation, which we refer the 
interested reader to for complete details. For our 
purposes, it will suffice to show the final result, 


which is Equation 7 in Abramo et al. ( 2007 1 : 


S, 


2 H - 


1 + Wi 


• 4 + 3 Wj Sj 

3 3(1 + Wj) 1 + Sj 


p k (l + w k )(l + 3w k )S k (l + S k ) = 0. (8) 


The subscripts j and k refer to different fluids in 
the system, e.g. matter and dark energy. In a 
scenario where the dark energy can clump, this 
becomes a system of equations. As we discuss be- 
low, we do not allow this possibility. Therefore, 
noting that w = 0 for matter, we get the follow- 
ing non-linear equation governing the behavior of 
matter perturbations: 


5 m + 2 HS m — — - — 4:TrGp m S m (l + 5 m ) — 0. 

3 1 + S m 

( 9 ) 

But, why exactly are we allowed to ignore the 
clumping in dark energy? 


The discussion about dark energy perturbations 
is often cast in terms of the effective speed of sound 
for dark energy, e.g., Bean & Dore (2004 1 . Typ- 
ically one may expect that for an adiabatic fluid 
with a constant equation state parameter w, the 
speed of sound is given by c 2 = 5P/5p = w. How- 
ever, when Wde < 0, clearly this becomes imagi- 
nary. suggesting a catastrophic instability, e.g., see 


Afshordi et al. (2005), and thus we must use a 


more general definition of the c. In other words, 
dark energy with constant w cannot be an adi- 
abatic fluid, implying w = ^ ^ |^. While the 
equation of state parameter can still give us infor- 
mation about the background evolution, the full 
action of the fluid is necessary to compute its effec- 
tive speed of sound: c 2 ^j = It turns out that 
for the simplest quintessence models c e ff = 1, al- 
though for more general actions c e ff could essen- 
tially take any value, e.g.,|Armendariz-Picon et al.| 
( 2001 ). 

As to the question of dark energy clumping, 
we know that pressure fluctuations propagate with 
the speed of sound c e ff. Therefore, dark energy 
should be smooth on scales smaller than its sound 
horizon ~ c e ff/H. As long as c e // ~ 1, all the 
collapsed structures at late times are much smaller 
than the sound horizon, implying that dark energy 
perturbations See should be negligible for their 
formation. 


3. Numerical Techniques 


Although Gunn & Gott (1972) showed that 
the perturbation equation with the assumption of 
spherical symmetry (also known as the spherical 
top-hat problem) can be solved analytically for 
the case of an Einstein-de Sitter (EdS) universe 
(A = 0), we are interested in cases where the 
cosmological constant/dark energy are non-zero. 
Therefore, a numerical solution is necessary. 

We built a code in C++ that utilizes the Runge- 
Kutta method for numerical solutions of ordinary 
differential equations. The full code can be found 


Prescod- Weinstein (2010). The code runs an 


instance of a loop for value of present hr the 
range of 0 to 0.7 (the currently measured value 
of dark energy density), which solves the second- 
order linear differential equation, Equation |A7[ as 
well as the full non-linear equation, Equation [6] 
In order to solve both equations, the solution to 
the differential equation for the background scale 
factor, Equation [4j is found. The curvature con- 
stant is calculated for the initial overdensity be- 
fore proceeding to a solution for both the linear 
and nonlinear differential equations. 


Collapse time is formally defined to be the time 
at which nonlinear overdensity, Snl, goes to infin- 
ity. However, we define collapse numerically by 
requiring a large value of nonlinear overdensity, 
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Snl = 200. The initial conditions are calibrated 
such that they provide the same results as the ana- 
lytical EdS model, namely that S^it-coiiapse) ~ 1.7. 
This is essentially done by assuming that at early 
times S scales linearly with the scale factor a, 
as ex pected for the growing solution to Equation 
( A7 ) in the matter era [^] 



Fig. 3.- This plot shows nonlinear vs. linear overdensi- 
ties in the Einstein-de Sitter Universe. The code produces, 
as expected, a Snl that diverges when Sl ~ 1.68. 


4. Results 

As stated in the Introduction, we wish to dis- 
cover the time (as a fraction of the collapse time) 
such that the variance of 5 Lit) is at a minimum, 
as we vary cosmologies. In doing so, we find the 
time in cosmological history where linear theory 
is most likely to accurately predict gravitational 
collapse, independent of cosmology. 

For simplicity, we chose to plot the relative dif- 
ference of 5 L{t) (in each cosmology) with respect 
to the fraction of time to the collapse time in the 
Einstein-de Sitter universe. However, it should 
be noted that the results are independent of this 
choice, and one could easily calibrate with respect 
to a universe with a non-zero A instead. Having 
found a common ground, all the data was searched 
for a single point in time (in units of collapse 
time) where the variance of sljBo) - 1 was at 
a minimum, where Q stands for different dark en- 
ergy models under consideration, whether a gen- 
eral fluid or a cosmological constant. 

The interpolations and variance computations 
were carried out using a script in MatLab, which 


, ~ 0.94, the 

t collapse 

apse) for different cosmologi- 


can be found in Prescod- Weinstein (2010). In the 
case of the simple cosmological constant models, 
Fig. 0 shows the result that at * 
variance in <5i(i/i co « 
cal constant models hits a minimum of 1.80 x 10 -9 . 
Values of Sl at range from 1.602 for Einstein-de 
Sitter 0 to 1.599 in the A = 0.7 universe. 

Computing the solution to Equation © re- 
quires some modifications to the code. Instead of 
scanning over different cosmological constant val- 
ues, this version of the code varies between con- 
stant values of w. Moreover, as current observa- 
tions, e.g., |Komatsu et ah] (2010) set a (very con- 
servative) upper limit of —1/3 for the value of w , 
we studied cases where w was smaller than —1/3. 

Because the cosmological constant is a special 
case of this scenario with w = —1, we are able 
to check the self-consistency of our methods (Eq. 
[6] vs. Eq. [9]) finding that both versions of the 
code produce the same results for a universe with 
Ha = 0.7 (approximately the universe that we live 
in). 



Fig. 4. — Relative change in SL(t/t co u apse ) in ACDM 
cosmologies = 0.1, 0.2, 0.7), with respect to the 

Einstein-de Sitter Universe. t co iia P se is calculated for 
spherical overdensities. The curves seem to intersect at 
t/tcoiiapse — 0.94, and a calculation of the point of mini- 
mum variance between the lines confirms this. 

Fig. § shows a similar comparison to that of 
Fig. 0, but with different values of equation of 
state, w , with £Ide — 0.7. Interestingly, we can 
again clearly see a point of minimum variance at 
t/tcoiiapse — 0.94, suggesting that this result might 
be quite independent of the dark energy model (at 
least within the spherical collapse approximation) . 


2 It is important to remember that <5 is a ratio of two numbers 
whose units are that of density. Therefore, 5 is dimension- 
less. 


3 This is below the predicted analytic value of 1.68 at collapse 
because while Sl = 1.68 is expected as Snl oo , we 

have set the collapse to occur at Snl — 200, resulting in a 
lowered collapse Sl- 
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Fig. 5. — Same as Fig. (J4ji. but with varying values 
of dark energy equation of state ( w = — 0.3, — 0.4, ..., — 1), 
at Vlf^(today) = 0.7. Again, curves seem to intersect at 
t / 1 collapse — 0.94, and, again, a calculation of the point 
of minimum variance between the lines confirms this. We 
begin at the value -0.3 because current data constrains the 
parameter to be smaller than this value. 


The variance at this minimum is 7.5 x 10 _ 7 . In the 
next section, we discuss the physical significance 
of this result. 


5. What does the Cluster Mass Function 
teach us about cosmology? 

At first, one might be puzzled by the fact that 
Sl happens to have almost the same value at 94% 
of the collapse time, independent of w or A, even 
though the linear approximation breaks down long 
before this point. In other words, why should non- 
linear collapse show such strong correlation with 
the linear evolution? This can be understood as 
the near cancelation of two different effects with 
opposite signs: 


1. With the exact same initial conditions, the 
presence of dark energy weakens the gravita- 
tional attraction near the turn around point, 
which in turn stretches the collapse time. 


2. The linear growth factor D(t ), which is the 
growing solution to Equation (A7), slows 
down as dark energy starts to dominate, 
since the Hubble friction 2 HS stops decaying 
( H -A const.), while matter density decays 
exponentially p oc a~ 3 oc exp(— 3 Ht). 


It turns out that for near ACDM cosmologies, 
these two effects nearly cancel each other, i.e. the 
linear growth is slowed down, but t co u apse is also 
longer, resulting in nearly the same value of Sl 
close to (i.e. at 94% of) the collapse time. 



log(l/cr) 


Fig. 6.— The expected time of minimum variance in 
linear overdensity, in units of ellipsoidal collapse time. Ob- 
serving the cluster number counts at t co i should tell us the 
linear overdensity of the collapsing region at £*, indepen- 
dent of the dark energy model. 


However, we should note that this result is 
specific to the spherical collapse scenario. On 
the other hand, real collapsing regions are far 
from spherical. Even though very rare peaks of 
a random gaussian field can be approximated as 
spheres, more abundant haloes could significantly 
deviate from sphericity, e.g., Bardeen et al. ( 19861. 


In addition, in a study of halo assembly histories, 


Angulo & White (2010) find that extending the 


PSF to account for ellipsoidal effects could sup- 
port a more accurate picture. For Einstein-de Sit- 


ter universe, Sheth et al. (2001) give a simple nu- 


merical fit for the impact of ellipticity on the linear 
collapse threshold, S ec : 


S ec — S sc [l + 0.31 x ct(M) 1 ' 23 ] 


( 10 ) 


where S sc ~ 1.686 is the spherical collapse thresh- 


old (but see Robertson et al. (20091 for the limita- 
tions of the ellipsoidal collapse formalism). Since 
S L oc t 2 / 3 in an Einstein-de Sitter universe, this 
implies that, for the same initial overdensity, the 
collapse time of an elliptical region is longer than 
that of a spherical region by a factor of: 


collapse, elliptical 
collapse, spherical 


~ [1 + 0.31 X cr(M) 1 ' 23 ] 372 . 

( 11 ) 
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In other words, we need to extrapolate the linear 
theory predictions in Figs. (4|5) farther beyond 
the point of intersection to actually hit gravita- 
tional collapse. Therefore, assuming that Equa- 
tion © is approximately independent of the dark 
energy model, it can be combined with the results 
of previous sections to show that the point of min- 
imum variance in 5l is shifted to smaller values of 

t/tcollapsei he.. 


0.94 


tcollapse [1 -)- 0.31 X cr(TFT) !-23J 3 / 2 


( 12 ) 


if we include the impact of ellipsoidal collapse. 
This result is shown in Fig. § , and demonstrates 
how measuring the mass function of clusters at a 
given era may tell us about the entire history of 
linear growth, and not just a snapshot at the time 
of observation (as is in the traditional universal 
mass function hypothesis). 

Inspired by the fitting formula Eq. §■ Equa- 


tion ( 12 1 leads us to propose a new universal clus- 


ter mass function: 


0.94 x t„ 


= „D [(1 + 031 x ^i. 23)3/2 

(13) 

where the actual mass function is related to 
f{cr',Q) through equation 0 , and t a is the time of 
observation at which D(t) is normalized to unity. 

In other words, Equation © suggests that the 
exponential cut-off in the cluster number counts 
at any time t is set by the linear density fluctu- 
ations <7* at an earlier time £*, which is set by 


Equation (12). As suggested by several recent nu- 


merical studies (see Introduction) f(a ; Q) depends 
on cosmology (denoted by Q), but we propose the 
functions g and h to be universal, while the de- 
pendence on cosmology (or dark energy models) 
only enters through D{t *). The explicit depen- 
dence of g and h on a at the time of observation is 
justified, as the value of a acts as a proxy for the 


asphericity of the collapsing region (Sheth et al. 


2001 ). 


Comparing to Eq. ([3]), we fix g and h as: 


h(a) = cD 2 


g(a) = A 
0.94 x t 0 


(?)' 


(1 + 0.31 X cr 1 - 23 ) 3 / 2 J q a=0 7 ’ 


(14) 


where ( A , a , b , c) = (0.186, 1.47, 2.57, 1.19) Q 


o.i 
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0 0.2 
log(l/<7) 
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Fig. 7. — The relative change in /(<r) at high redshifts, 
compared to z = 0. The colored regions show the simu- 
lated results from|Tinker et al. ( |2008[ ). Curves in top panel 
show our analytic prediction without ellipticity corrections, 
while curves in the bottom panel include the ellipticity cor- 
rections (Eq. |13[ ). The solid, dotted, short-dashed, and 
long-dashed curves refer to z = 0, 0.5, 1.25, and 2.5 respec- 
tively. 

Fig. 0 compares the evolution in the func- 


tion /(cr; Q ) in N-body simulations of Tinker et al. 


4 We recognize that these choices for g and h are not unique, 
but as Fig. 0 demonstrates, these are sufficient to fit 
current simulations 


(2008), with our prediction from Eqs. (IsjM). We 
see that while the result from spherical collapse, 
f* = 0.94 x t co u pase , underpredicts the evolution 
(top panel), Eq. (13), which includes the ellip- 
soidal correction to the collapse time (Eq. [12) ), 
can successfully reproduce the evolution in /(er; Q ) 
(bottom panel). Notice that, while we use the 
same er-dependence as Tinker et al. ( 2008 1 at z = 0 
to fix functions h(a) and <?(< r), our predicted red- 
shift evolution is based on Eq. (13). Therefore, 
comparison with their simulated evolution pro- 
vides an independent test of our proposal. 

In the end, we should point out a subtlety in 
comparing our predictions with simulated cluster 
mass functions. The most widely used algorithm 
to identify haloes in N-body simulations, known as 
friends-of- friends (FoF), is based on linking neigh- 
boring particles. However, this method has the 
tendency to bridge between two close-by haloes, 
and thus counting them as one. An alternative 
method, which is known as the Spherical Over- 
density method, (SO) identifies haloes by search- 
ing for spheres of large overdensity (200 or more). 
The SO haloes are more closely related to the ob- 
servational notion of clusters, as large overdense 
compact regions of the Universe. Consequently, 
we expect our predictions based on spherical (or 
ellipsoidal) collapse to be relevant for SO cluster 


mass functions, such as in Tinker et al. (2008). 


Indeed, the deviations from universality show op- 
posite trend with mass for FoF clusters, with more 
(less) massive clusters having closer to (farther 


from) universal behavior ( 

Reed et al.|2007j Crocce 

et al. 2010 

Bhattacharya et al. 

2010). We note 


that predictions from FoF and SO methods are 
divergent, leading to substantive challenges to our 
ability to compare models with cluster observa- 
tions. Further discussion of this topic, however, is 
beyond the scope of this article. 


6. Conclusions and Future Prospects 

We have presented a study of non-linear grav- 
itational collapse in different models of dark en- 
ergy/cosmic acceleration. In particular, we crit- 
ically examined the correlation between the lin- 
ear growth of fluctuations and the emergence and 
statistics of collapsed objects (such as dark mat- 
ter haloes or galaxy clusters). First, we focused 
on the collapse of spherical overdensities, and dis- 


covered that they all have the same linear ovder- 
density (~ 1.50), at ~ 94% of the time of col- 
lapse/virialization, independent of the density or 
equation of state of dark energy. We then used a 
simple prescription to include the impact of el- 
lipsoidal collapse in our finding, and then used 
this result to propose a new universal mass func- 
tion for galaxy clusters/dark matter haloes. Our 
semi-analytic predictions are consistent with the 
observed evolution in mass function of haloes in 
N-body simulations. 


Future work will include the adaptation of this 
prescription to study models with dynamical equa- 
tions of state, such as quintessence or modified 
gravity models. A particularly challenging appli- 
cation would be to the gravitational aether/black 
hole model ( Prescod- Weinstein et al. 2009). Be- 
cause of the way dark energy is sourced in the 
gravitational aether model, there are unique tech- 
nical challenges associated with properly describ- 
ing structure formation in its presence. Early 
Dark Energy (Grossi & Springel 2009) cosmolo- 
gies also provide an interesting example of a 
quintessence model where dark energy signa- 
tures in structure formation would be particularly 
strong.) 

Finally, we recognize that a more systematic ap- 
proach to the question of universal mass functions 
should be possible, and given the level of theoret- 
ical and observational activities in this field, it is 
unlikely that the present work provides the last 
word on this subject. For example, the result will 
be sensitive to the shape of the primordial power 
spectrum (Bagla et al.|2009), and can be extended 


to the conditional mass function (e.g., Neistein 


et al. t2010 ). However, the novel (albeit trivial) ob- 


servation of this work is that measuring the clus- 
ter mass function will teach us about the entire 
history of linear growth. This is in contrast to 
many previous cosmological applications of cluster 
mass functions, which assume a one-to-one corre- 
spondence with the linear a{M) at the time of 
observation. Similar to a multi-level archeologi- 
cal excavation, dark matter haloes can be thought 
of as artifacts of linear growth. As Fig. § and 
Eq. (13) suggest, the number of smaller haloes 
(with larger u(M)) can teach us about the early 
evolution of linear growth, while the bigger haloes 
(with smaller er(M)) tell us about its more recent 
history. 
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We thus speculate that this perspective can 
eventually lead to yet untapped information about 
the nature and history of cosmic acceleration, es- 
pecially as the releases of several large scale clus- 
ter surveys such as Atacama Cosmology Tele- 
scope (ACT) ( Menanteau et al.|2010 ), South Pole 
Telescope (SPT) (Andersson et al. 2010), Planck 
(Geisbiisch & Hobson |2007 ), and Red Sequence 
Cluster Sequence 2 (RCS2) (Yee et al. 2007) are 
now looming on the horizon. 
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A. Linear Perturbations 


The linear perturbation theory that is used to describe structure formation can be derived via a fluid 
picture. We use a Newtonian treatment because when density perturbations are small, the gravitational 
potential will be nonrelativistic ( Peebles(l 993 ). The standard equations of fluid dynamics are reviewed here. 
First, the continuity equation: 

r) r 

(Al) 


f + v. W = o, 


where p is the matter density and v is the fluid velocity. The Euler equation is: 

dv VP 

— + (vV)v= --Vip, 

where <fi is the gravitational potential. Finally, the Poisson equation is: 

V 2 0 = AttGnP- 


(A2) 


(A3) 


Introducing the comoving coordinates x = , we can write the velocity in terms of the comoving coordinates 

and the scale factor a(t): 

rlr . n r 

(A4) 


- dr . _ •_ a _ _ . r . 

v = — = ax + xa = —r + u( — ,t) 
at a a 


where u = xa(t) is the peculiar velocity. Note that at constant r , ^ = 0 so ax + ax = 0, giving us ^x = — 
This leads us to the following relation: 


dt ’ 


( 9p(r, t ) 

l dt 


= ( dp{r, t) 

r V dt 

dp{r, t) 
dt 


dx dp 
dt dx 


- -x- V x p{x,t). 


(A5) 


We now wish to rewrite Equation |A1| in terms of the comoving coordinates, which essentially means 
replacing the partial differential with the modified one from Equation |A5| Moreover, what we are really 
interested in is the development of relative deviations from the mean density, or the density contrast: <5 = 
| — 1. Thus, we write p = p( 1 + S) and assuming that p is the mean density of regular matter, we expect 
p oc a~ 3 . This gives: 

°=(f)pi v ' (1+i)c (A6) 

We make similar transformations for the Poisson and Euler equations. Next, we drop higher order terms (e.g. 
0(u 2 ) or u6). We differentiate the linearized continuity equation and take the divergence of the linearized 
Euler equation. This gives us a differential equation that depends entirely on S and not on u: 

5 + 2—8 = AttGnpS. (A7) 

a 

The linear growth factor D(t) is defined as the growing solution for 5 in this equation. 
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